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Abstract. The thermochemistry of H 2 and HD in non-collapsed, non-reionized primordial gas up to the end of the 
dark age is investigated with recent radiation-matter and chemical reaction rates taking into account the efficient 
coolant HD, and the possibility of a gas-solid phase transition of H 2 . In the standard big-bang model we find that 
these molecules can freeze out and lead to the growth of flakes of solid molecular hydrogen at redshifts z ~ 6 — 12 
in the unperturbed medium and under-dense regions. While this freezing caused by the mere adiabatic cooling 
of the expanding matter is less likely to occur in collapsed regions due to their higher than radiation background 
temperature, on the other hand the super-adiabatic expansion in voids strongly favors it. Later reionization (at 
z fs 5 — 6) eventually destroys all these H 2 flakes. The possible occurrence of H 2 flakes is important for the degree 
of coupling between matter and radiation, as well as for the existence of a gas-grain chemistry at the end of the 
dark age. 
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1. Introduction 

The period of time between the last scattering of the back¬ 
ground radiation photons and the formation of the first 
bound objects (i.e., the “dark age”) still presents inter¬ 
esting unknowns even about the main processes that hap¬ 
pened then, as described below. 

The studies of the physical and chemical conditions in 
the post-recombination Universe clearly established that 
trace amounts of molecules such as H 2 , HD and LiH 


were present in the primordial gas (Lepp & Shull 1984; 
Dalgarno & Lepp 1987; Latter & Black 1991; Puy et 
al. 1993; Stanch et al. 1996, 1998; Galli & Paha 19981 ). 
A recent and comprehensive revie w of primordial chem¬ 
istry is given by Lepp et al. ( 20021 ). 

Because of its rotational levels, H 2 is an important 
coolant of interstellar gas from 10^ to 10^ K. Despite its 
lower abundance, the molecule HD is also important, since 
its allowed dipole rotational transitions provide strong 
channels coupling radiation and matter at temperatures 
below 10^ K, well below that of the lowest H 2 quadrupole 
transitions (at 512 K). Thus primordial molecules can play 
an important role in the post-recombination Universe, 
already in the uniform, non-collapsed fraction of the 
Universe. 


Molecules can play an even more crucial role in the 
first collapsing structures appearing at temperatures typ¬ 
ically below a few hundred K, since they offer the possi¬ 
bility to efficiently cool gaseous proto-clouds (see Lahav 


1986; Puy & Signore 1996, 1997; Uehara et al. 1996, 2000 


Abel et al. 1997, 2000). However a definitive understand¬ 
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ing of the collapse of the first structures is still lacking, 
because the coupling of gravity, chemistry, and radiation 
constitutes a formidable non-linear system for which not 
all equations are presently well known. Since the system of 
equations is “stiff”, i.e., exponentially sensitive to pertur¬ 
bations, changes in the conclusions can eventually arise if a 
s ingle factor is modified. For instance Combes & Pfenniger 
( 1998 ) show that the fraction of molecular hydrogen can 
reach unity if the collapsing structures grow not as smooth 
isolated spheres, but similarly to the fractal structures ob¬ 
served in the galactic ISM. 

Zwicky ( 1959 ), Takeda et al. (|l969 ), Hirasawa et 
al. ( 1969 ) and Matsuda et al. ( 1969| ) pointed out early that 
molecules could appear during the post-recombination pe¬ 
riod. The mechanism by which molecules form in the post- 
neutralized Universe is traditionally assumed to be differ¬ 
ent from the chemistry of the interstellar medium, because 
the primordial chemistry is supposed to remain a pure gas- 
phase chemistry. 

The first natural extension of molecular physics toward 
larger matter aggregates is the possibility of a gas-solid 
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phase transition by molecular hydrogen in the primordial 
gasQ. This is a central question, because the occurrence 
of H 2 ice grains, or more likely fluffy structures such as 
flakes, may radically change the chemistry via grain sur¬ 
face reactions, as well as the matter-radiation coupling via 
photon-grain interactions. 


Solid molecular hydrogen^, the simplest molecular 
crystal, is rather special among the molecular crystals be¬ 
cause at low pressure the weakly anisotropic H 2 molecules 
are free to rotate and vibrate in the lattice. Further, as 
for helium, due to its large (comparable to the binding 
energy) zero-energy motion, quantum effects are impor¬ 
tant; such crystals are known as quantum crystals. Solid 
hydrogen has been studied in the laboratory for at least 70 
years, and used in several cryogenic and industrial applica¬ 
tions, so many of its properties are today well established. 
A comprehensive review of solid molecular hydrogen is 
given by Silvera ( 1980 ). 


In this paper only first order effects associated with 
the condensed phase of molecular hydrogen are consid¬ 
ered. Much more work will be required to describe other 
possible effects occurring in the non-linear regime follow¬ 
ing gravitational instability. As it turns out, much depends 
on the temperature of the matter, and whether matter can 
be cooled in the interval ~ 1 — 2 K by mere adiabatic ex¬ 
pansion. 


Mechanical cooling is usually discarded because the 
overwhelming number of photons in the cosmic back¬ 
ground radiation is thought to couple well enough matter 
and radiation with the few residual electrons by Thomson 
scattering (see, e.g., the clear discussion by Longair 1995. 
p. 428). Here we calculate in detail the thermo-chemical 
evolution with the most recent published rates when the 
matter temperature cools as ^ R~^, faster than the ra¬ 
diation temperature ^ R~^ due to the universal expan¬ 
sion {R is the Universe scale factor). The adiabatic cool¬ 
ing ~ is found to be effective below redshifts ^ 300, 
leaving the possibility to cool non-collapsed, non-reionized 
matter to very cold temperatures, until most of intergalac- 
tic matter is eventually reionized. Current models set the 
almost complete reionization around redshifts 2^5 — 6 
(Razoumov et al. 2002). 

The outline of this paper is as follows. In Sect. 2, the 
method of our calculations and the chemistry is described. 
In Sect. 3, results about the revised thermochemistry are 
presented, including the calculation of the phase transi¬ 
tion to solid H 2 . The implications are discussed briefly in 
Sect. 4. Finally, our work is summarized in Sect. 5. 


^ No corresponding gas-solid transition exists. Below 
5.2 K, the critical point of “^He, down to 0 K the “^He gas-liquid 
phase transition occurs at much higher pressure than for the 
H 2 gas-solid transition. 

^ Not to be confused with metallic hydrogen, a high pressure 
form of hydrogen, found for example inside Jupiter. 


2. Model 


In most interstellar environments, association catalyzed 
on grain surfaces is thought to dominate the formation 
of H 2 . Thus, molecular hydrogen is generally assumed 
to form by hydrogen atom recombination on cold grains 
in dense molecular clouds (Hollenbach & Salpeter 


1970 


1971 ). At the epoch of recombination where a total ab¬ 
sence of dust grains appears justified, surface reactions 
are ignored. The two standard and principal gas-phase 
processes are characterized by 1) the process: 


H + H+ 
H++H 


H+ 


hv, 

H2, 


and 2) the H process: 


H-he" 
H" +H 


H -I- hv, 
e- +H2. 


( 1 ) 


( 2 ) 


However, these reactions are coupled with other reactions 
because the main chemical species available from primor¬ 
dial nucleosynthesis are H, H+, D, D+, He and Li. In par¬ 
ticular this simple chemical network must be coupled with 
the reactions which lead to the formation of HD molecules, 
because HD is an efficient coolant below ^ 100 K. 


Recently Galli & Palla (2002) gave an updated re¬ 
view of the deuterium chemistry of the post-recombination 
Universe from their exhaustive work on the primordial 
chemistry (Galli & Palla 199^ ). The formation of HD in 
the primordial gas follows the two routes: 


D+ -bH2 
HD+ +H 


H+ 

H+ 


HD, 

HD. 


The radiative association 


H-bD 


HD -b hi^. 


( 3 ) 


( 4 ) 


is very slow in our context. 

The chemical reaction rates depend on the tempera¬ 
ture of the species or of radiation (for the photoprocesses), 
and of the density. Thus the chemical kinematics must be 
coupled with the hydrodynamic equation in the framework 
of a cosmological model (here the Friedmann model). 

A rigorous approach would be to consider four fluids: 
neutrals, positive ions, electrons and photons with their 
corresponding temperatures in the hydrodynamic equa¬ 
tions. However Flower & Pineau des Forets ( pOOO ) showed 
a very small difference between the temperature profiles 
of the neutrals, ions and electrons in the primordial gas. 
In fact, the coupling through Goulomb scattering between 
the electrons and the ions and the neutrals is sufficiently 
strong at the post-recombination epoch. In the following 
Tin, is the common temperature of the neutrals, ions, and 
electrons, which we call matter temperature. 

The temperature of the cosmic background radiation 
(GBR) is given by: 


dT, 

dt 


= -H{z)T,, 


( 5 ) 
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Ho 

67 Km s ^ Mpc ^ 


0.24 

^r,0 

9.265 X 10"® 

[D] 

3.3 X 10“® 

f^m,0 

~ 0.29990 

[Li] 

2.1 X 10"“ 

Ub,o 

0.0535 

[H2] 

10-40 

Uk,o 

0 (flat Universe) 

[HD] 

10-50 

Ua,o 

0.7 

[LiH] 

10-50 


Table 1. Standard cosmological model (S-model): parameters and initial relative abundances at 2 : = 10'*. 


a relation which has been observationally confirmed by 
Srianand et al. ( ^000 ), at the redshift z = 2.138, with 
absorption lines from the first and second fine-structure 
levels of neutral carbon atoms in an isolated cloud of gas. 

H{z) is the Hubble parameter at the redshift z, defined 
from the scale factor R, and Hq (the Hubble constant) is 
its present value: 

H{z) = R/R = (6) 

1 

Hq r2i.,o(l + .z)^ + + + + ^^A,0 j 

where the are constants with the present values, con¬ 
strained by the Friedman-Lemaitre condition: 


f ^ r ,0 + f ^ m ,0 + ^^^,0 + f ^ A ,0 — 1 - 

the radiation density parameter, is given by: 


(7) 




dt 


— —‘2H{z) Tjji -|- dlCompton + ^mol- 


(9) 


4- 


Compton 


8 axa 

3 TOpC 


T^{T,-T^)x,. 


( 10 ) 


function and gives a slight contribution to the evolution 


of matter temperature (see Puy et al. 1993). 

This series of differential equations is completed by the 
redshift equation: 


dz 

dt 


= -H{z){l + z), 


( 11 ) 


and by the equation of the numerical density of baryons 
n: 


dn 

dt 


= —3H{z) n, 


( 12 ) 


which just describes the dilution of the conserved number 
of baryons in the expanding volume. 

Each chemical species rig has a similar behavior, com¬ 
pleted with the chemical network: 


Q Tjf ^ 

— = -m{z)n, 


/dns\ 

V dt / chem 

where the chemical network 


where a = Aiojc is the radiation constant related to 
the Stefan-Boltzmann constant cr, Tsk = 2.726 K is the 


(13) 


has the form: 


chem 


present radiation temperature (Mather et al. 1994 ). The 
neutrino contribution to the radiation density for three 
massless neutrino types is noted /i,, G is the gravitational 
constant, and c is the speed of light. is the matter 

density parameter (i.e., baryons and other non-baryonic 
matter), is the curvature parameter, and Ha.o is 

the mass-energy density parameter associated with the 
Einstein’s cosmological constant A, or with quintessence 


( drig 


It.) — ^ OlijTLiTtj PijTls- 

V dt /chem ■> J 


(14) 


(Wetterich 198S, Peebles & Ratra 2002). In our context, 
Ha.o will be held constant., fixed for z=0. 

The evolution of the gas te mpera ture Tm is governed 
by the equation (see Puy et al. 1993|): 


4^ Compton characterizes the energy transfer from radiation 
to matter via Compton scattering of CBR photons on elec¬ 
trons (Peebles 1968): 


<7t defines the Thomson cross-section and rrie the electron 
mass. Xe is the electron abundance (i.e., the ionization 
fraction), 'kmoi is the energy transfer via excitation and de¬ 
excitation of molecular transition. However, in the post¬ 
recombination context this thermal function is a heating 


The first terms of the sum define the s formation process 
(* + J ^ s) associated with the reaction rate , the 
second term characterizes the destruction process (s ^ 
i+ j) with the reaction rate Ptj. Some reactions rates 
and (3ij are given in the Appendix. 

2.1. Initial conditions 

With stiff systems of equations, it is obviously important 
to use the most accurate available physical constants and 
parameters. In the context of the primordial chemistry the 
knowledge of the abundances of the primordial nuclei is 
crucial. Prom the latest determination of primordial he¬ 
lium (i.e., Yp ^ 0.24) and deuterium (D/H ^ 3.3 x 10“^) 
abundances given by Tytler et al. ( pOOOD we can deduce the 
baryon density, which is in good agreement with the cos¬ 
mic microwave background meas ureme nts of Boo meran g 
and Maxima (de Bernardis et al. 200C, Jaffe et al. 2000 ): 
flb,o ~ 0.024. These experiments find a peak in the an¬ 
gular power spectrum of the microwave background, con¬ 
sistent with that expected for cold dark matter models in 
a flat Universe (i.e., flK,o = 0)- 

The observations of type la supernovae (Perlmutter 


et al. 199£, Riess et al. 1998) suggest that the Universe 
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may be presently dominated by an additional dark en¬ 
ergy. Combined observations of type la supernova e, cos - 
mic microwave background anisotropy (Jaffe et al. 2000) 
and cluster evolution (Bahcall & Fan |1998D for which the 
results have been done in the form of likelihood contours 
in the flm.o and Oa.o plane, yield 


f^m.o ~ 0.3 and flA,o ~ 0.7. 


(15) 


1-5 


In our case flm.o ~ 0.2999 and flr.o ~ 9.265 x 10” 
order to verify the Friedman-Lemaitre condition (Eq. 

We adopt the Hubb le con stant Hq = 67 Km s“^ Mpc“^ 
given by Freedman ( ^000 ). We summarize the parameters 
of our standard model (S-model) in Table 1. 

In other models (see discussion) we will vary only the 
value of the baryonic fraction flb.o- Cosmochemistry evolu¬ 
tion is not sensitive to the Ha.q parameter. Molecular for¬ 
mation appears in the redshift range [100,1000] whereas 
Ha.o plays a role in the Hubble parameter, see Eq. Jlat 
redshift below z ~ 5 (see Carroll, Press & Turner 1992|) . 
We will fix this last parameter to the value Ha.o = 0.7. 


2.2. Chemical evolution 

The numerical integration of the coupled chemical and hy- 
drodynamical equations is an initial value problem for stiff 
differential equations. We use the ch emical network devel¬ 
oped by Galli & Palla ( 1998 , 2002 ), except the hydrogen 
and deuterium recombination which is calculated from the 


reaction rates given by Abel et al. ( |1997|) , see Appendix. 

We start our integration at redshift Zi = 10^. At this 
value it is particularly important to take into account^ the 
radiation contribution Hj. o in order to calculate the age of 
the Universe t, at redshift z, : 


— 1/2 fjrp 

— (1 
X 


I [°° \ 

ti — I n," qX -t-f2A,0 

^0 Ji+z2 

The time U gives a timescale for the initial integration of 
the se t of E qs. d-©. which is solved by Gear’s method 
(Gear |l97l|) . 

In Pig. 1(a) we show the helium and hydrogen recombi¬ 
nation which give the abundance of the free electrons. The 
slight decreasing of the electronic fraction Xe at z ~ 2485 
corresponds to the helium recombination. At the redshift 
z ^ 1360 the medium becomes neutral (hydrogen recom¬ 
bination: neutral hydrogen is dominant). In all of these 
cosmological models the fraction of electrons, Xe, is low 
at z = 1. For the S-model Xe = 3.04 x 10~ ^ is i n close 
accordance with the results of Galli & Palla ( 1998| ) where 
Xe ~ 3.02 X 10“^. Before and around recombination, the 
thermal history depends on the tight coupling between 
radiation and matter resulting from Thomson scattering. 
Well after recombination, the free electrons quasi vanished 


® In our cosmological model, the equilibrium redshift (i.e., 
between matter and radiation density) is U-Zeq = IIm,o/Dr,o = 
-Dm,o ~ 3235. This value justifies that it is nec- 


3gg 

SttG aT* 


1 

1 + /..' 


essary to consider the radiative component. 





Fig. 1. a): Helium and hydrogen recombination for a fiat 
Universe with Hq — 67kms“^Mpc“^ and Ua,o = 0-7. 
The S-model corresponds to the baryonic fraction klb.o = 
0.0535, the M-model to Ub o = 0.1, and the B-model to 
nb,o = 1. b): Thermal decoupling between radiation and 
matter in the S-model. c): Gomparison of the H 2 and 
HD density in the S-model. After a transient growth H 2 
and HD follow the evolution of the neutral density, which 
means that the fractional abundance becomes almost con¬ 
stant below z ~ 500. 
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-i.e Xe ~ few [10“^ — 10“®], see Fig. 1(a). The evolu¬ 
tion of the mean temperature of radiation and of matter 
is shown in Fig. 1(b). The subsequent evolution of mat¬ 
ter temperature is mainly caused by the dilatation of the 
Universe (first term in the second member of Eq. de¬ 
spite the very high numerical photon density at low red- 
shifts (see Longair 1995, Partridge 1995). We notice that 
the thermal decoupling is effective at redshift Zdec ~ 498 
(we define Zdec the redshift below which the ratio Tm/Tr is 
lower than 0.99). The non-collapsed matter temperature 
can drop much below the radiation temperature, unless 
some other coupling channels not included in our model 
are effective. The possibility of cooling matter well below 
radiation is crucial for the possibility to form H 2 flakes. 
In Fig. (Ic) we have plotted the evolution of the density. 
After the thermal decoupling the total density is mainly 
due to neutral hydrogen. The H 2 density rises through two 
processes: the first step corresponds to the FtJ channel and 
the second one to the H“ channel. The HD density rises 
also in two steps with a transition around z ~ 100, the 
(D+, H 2 ) channel and the (HD+, H) channel. 


3. Thermochemistry 


Fifty years ago, van de Hulst ( 1949 ) mentioned the 
possibility of the existence of solid H 2 in the inter¬ 
stellar medium, which implies the possibility of a gas- 
solid phase transition. This possi bility was elaborated 
by Wick ramasinghe & Reddish (|196^ ) and Hoyle et 
al. ( 1968 ). In the context o f bary onic dark matter as cold 
H 2 Pfenniger & Combes (1994) and Wardle & Walker 
(1999) found also that solid H 2 can play an important 
role. 

This form of molecular hydrogen was already well 
known in the laboratory for a long time. For example, 
Harold Urey in his Nobel prize lecture (1935) described 
how this phase transition was used to easily separate the 
H 2 and D 2 isotopes, thank to their slightly different triple 
points. 

From cryogenic industry tables (L’Air Liquide 1976D , 
the sublimation curve of para-hydrogen can be well fit 
with the following equation (Pfenniger & Combes 1994): 


Psat 


= 5.7x 10 


)”ry^exp(]-5^]) 


K cm 


-3 


in the range of matter temperature Tm between 1 K and 
the triple point te mper ature at T 3 = 13.8 K. According 
to L’Air Liquide ( 1976 ) p. 890, below about 20K over 
99% of normal H 2 is in the para form. The asymptotic 
mixture containing the maximum of ortho-H 2 takes place 
between about 20 to 300 K, above which the ortho/para 
ratio converges toward 3. However, Flower & Pineau des 


Forets (200C), investigating the ortho/para H 2 ratio, found 
that this ratio decreases to a frozen value of about 0.25 for 
z < 20 due to the low density of the main conversion agent, 
H+. We will see such a different ortho/para ratio has a 
very small effect in the result about the H 2 condensation 
redshift. 




Fig. 2. (a): Comparison of the ratio between the gas 

pressure and the saturation pressure Psat, for H 2 -HD gas 
between IK and T 3 = 13.8K (i.e., range of redshift be¬ 
tween 1 and 30), Zflake is the redshift below which flakes 
form; B-model: Ub,o ~ 1, M-model: Ub,o ~ 0.1, S-model: 
f^b.o 0.0535 (for all of these models Ua,o = 0.7). For 
each model, the solid curves corresponds to molecules with 
allowed rotation, while the dashed curves to non-rotating 
molecules, b): Comparison of the saturation pressure Psat 
with the H 2 partial pressure Ph 2 j for fh® S-model as a 
function of temperature. 


If we suppose an ideal (H 2 -HD) gas the ratio ^ 
between the H 2 -gas pressure and the saturation pressure 
becomes: 

^ = Elk ~ 2.63 X 10“^^nH2T“3/2exp[^i(^Y (18) 

where riHa is the H 2 density (in cm“^, see Fig. 1(c)). 
Although HD has a dipole moment when H 2 has 
quadrupolar transitions, the thermodynamics of HD gas 

^ In this Sect., ph 2 = |fiH 2 7m where we neglect the rotation 
and vibration of the molecule, a valid approximation in the H 2 
solidification range 0 — 13.8 K. 





























akes in the early Universe 



Fig. 3. Evolution of temperatures of matter and radiation. 
The critical temperature is reached at redshift 6; below 
this redshift H 2 could stick to the surface of ice. 


is assumed to be very close to H 2 (without distinction 
between H 2 and HD). 

In Fig. 2(a) the ratio ^ is plotted for three different cos¬ 
mological models: 1) S-model: the standard cosmological 
model, 2) M-model: model with flb.o = 0-1) 3) B-model: 
model with flb.o = f^m,o ~ 0-3- For all of these models 
= 0.7). We define znake as the redshift for which be¬ 
low this boundary the ratio ^ is greater than 1, and thus 
leads to the possibility to freeze the (H 2 -HD) gas. 

We find that Zflake ~ 12 for the B-model, Zflake ~ 7 for 
the M-model, and Zflake ~ 6 for the S-model. The redshift 
-^flake is larger in the B-model because the intrinsic baryon 
density is more important. 

In Fig. 2 we show also the very slight influence of the 
molecule rotation factor, 5/2 or 3/2, in the H 2 pressure. 
The redshift Zflake is not essentially affected. The same 
remark is valid when the ortho/para ratio is taken as 0.25 
instead of 0, since ph 2 is lowered by only a factor 0.75. 
In Fig. 2(b) we compare the pressure Psat and the partial 
pressure ph 2 (S-model) below the triple point. Obviously, 
the insensitivity of Zflake to variations of ph 2 results from 
the steepness of the H 2 sublimation curve. 


4. Discussion 

What we have shown above is that below a critical redshift 
a necessary condition to form H 2 solid grains is fulfilled. 
However, this condition is not sufficient, in reality several 
complications may prevent or favor the formation of solid 
H 2 . We discuss below some of these complications, as well 
as some consequences in case solid H 2 does form. 


4.1. Coupling of chemistry and H 2 flakes 


Sandford & Allamandola (1993) showed experimentally 
that H 2 molecules can stick to the surfaces of usual ices, 
and argued that gas-phase H 2 freezes out onto classical 


dust grains. Thus H 2 -containing ices may be fairly com¬ 
mon in dense molecular clouds. 

On the other hand, solid hydrogen is a quantum solid 
with a large zero-point energy (Silvera 1980), which means 
that contrary to usual dust grains, atoms or molecules 
adsorbed on the surface of H 2 ice do not stick at fixed 
lattice sites but move easily around to go to the lowest 
free energy configuration. The consequence is that surface 
chemical reactions may be thus greatly accelerated on H 2 
ice. 

Atoms or molecules stick on the H 2 ice surface below 
a critical temperature which is estimated by Sandford & 
Allamandola: 


Tc.r\t — 


100 


56.46 — ln(nH2) 


K. 


(19) 


In Fig. 3, we compare the evolution of the radiation 
and matter temperatures with the evolution of the critical 
temperature. We see that the critical redshift ZcHt = 6j ke., 
also the critical redshift for which the temperature of (H 2 - 
HD) gas is below the critical temperature (Eq. The 
value of Zcrit is remarkably close to the threshold redshift 
leading to the formation of H 2 flakes (Figs. 2(a) and 2(b)), 
which is not unexpected since condensation and sticking 
phenomena are both closely related to the inter-molecular 
force. 

Thus, we can imagine that once the first flakes form, 
the formation of H 2 and HD may be enhanced by neutral 
hydrogen reactions on the flake, which in turn increase the 
H 2 vapor pressure, increasing the flake growth. 

Sandford & Allamandola ( 1993 ) speculated on the 
possibility that pure H 2 ices might exist in intergalactic 
and extragalactic space. This last hypot hesis had received 
some atten tion previously (Lee et al. 1971 or Hegyi & 
Olive 119861) . 


4.2. Coupling of background radiation and H 2 flakes 

However, in the still neutral and not collapsed regions of 
the Universe, the background radiation is warmer than the 
adiabatically expanding matter. If H 2 flakes form, a new 
channel opens for coupling matter and radiation, since 
solids possess numerous phonon bands able to interact 
with infrared radiation. 

The most prominent absorption bands in solid H 2 oc¬ 
curs in the interval 1.89 — 2.50 pm (e.g.. Allin et al. 1955 ), 
and at higher frequency in the interval 1.06 — 1.24 /im 
(Varghese et al. 1987 ). They are in part caused by the 
vibrational and rotational quadrupolar absorption as for 
the isolated molecule, and by the main vibrational mode 
between two adjacent molecules in the lattice. 

A particularity of solid H 2 is that the molecules in the 
crystal at low pressure are almost free to vibrate and ro¬ 
tate, molecular transitions of low quantum numbers being 
allowed. Therefore it may be useful to mention a few other 
transitions for free molecules in the far infrared. Zwicky 
(1959) mentioned a line at 85 pm associated to the H 2 
ortho-para transition, but this transition is very slow, of 

















Daniel Pfenniger and Denis Puy: H 2 flakes in the early Universe 


7 


order of 300 yr. Particularly interesting for the dark age 
context is the HD dipolar rotational transition around 
114/rm (e.g., Lee et al. 1988) which, despite a low HD 
abundance, is relevant since the transition is dipolar, and 
the background radiation particularly intense in the far 
infrared. 


Thus any coupling of radiation with solid H 2 should 
increase the matter temperature, which acts to sublimate 
H 2 flakes, a self-destructive process. At any redshift an 
equilibrium temperature should be reached where the con¬ 
densation and sublimation rates are equal. Several uncer¬ 
tainties prevent us from estimating these time-dependent 
rates, such as the radiation-flake coupling which depends 
on the size, geometry and crystalline structure of the H 2 
flakes, and also the chemistry coupling mentioned above. 
More “cosmological” laboratory experiments simulating 
aspects of the physico-chemical conditions during the dark 
age would be useful. 


4.3. Heating from latent heat 

The formation of H 2 ice liberates molecular binding en¬ 
ergy in the form of heat, the latent heat. In the range 
1 — 3 K the latent heat amounts to 92 — 99 kcal kg“^ (L’Air 
Liquids 1976), or 22 — 62 times the molecule kinetic en¬ 
ergy f/cT. Since the H 2 abundance is then of the order of 
10 “®, clearly the latent heat can easily be shared with the 
neutral gas without causing any significant temperature 
increase. 


4.4. Destruction of H 2 flakes by reionization 


At redshifts 5 — 9, one expects that a growing fraction of 
the non-collapsed Universe is filled by the UV radiation 
produced by the first collapsed and light emitting struc¬ 
tures, stars and quasars (see, e.g., the models of Razoumov 
et al. 2002). This means that at redshift 6—7 possibly most 


of the diffuse gas has already passed through a ionization 
front, and is reheated to high temperatures. So the forma¬ 
tion of flakes may be partly or totally suppressed if the 
reionization occurs early enough, and propagates in ev¬ 
ery corner of the Universe. Voids are the last regions to 
be reionized, for reasons well discussed by Razoumov et 
al. (| 200 ^ ): in brief, UV radiation percolates through ion¬ 
ization fronts and thus ionizes the dense regions first. The 
number of UV photons necessary to reach and ionize voids 
is therefore much larger than, as earlier assumed, for a UV 
optically thin medium bathed by a general radiation field. 
The result is to delay the full reionization of the Universe, 
leaving some more time for H 2 freezing, particularly in 
voids. 


4.5. Formation of H 2 flakes in voids 

An additional mechanism enhancing H 2 freezing occurs in 
regions expanding faster than the Universe, i.e., in voids, 
since there the adiabatic cooling is even stronger. The es¬ 


sential factor favoring H 2 freezing is the temperature. If 
the void matter reaches a temperature around 1 — 2 K the 
formation of H 2 flakes is very much facilitated because the 
saturation pressure drops extremely fast, by over 27 dex 
in this range, as can be derived from Eq. |^, or Fig. 2(b). 
Since in an adiabatic mono-atomic gas T oc a region 
that is 10 times underdense with respect to the average 
would be, for 7 = 5/3, 10^/^ = 4.6 times colder. Since 
T„i oc (1 -|- z)^, the redshift at which H 2 reaches the criti¬ 
cal freezing temperature would about double, bringing it 
at a much darker epoch z > 12 . 


4.6. Formation of H 2 flakes in collapsing structures 


On the other hand, if gas radiative cooling oc rF is suf¬ 
ficiently strong, in collapsing regions eventually the tem¬ 
perature increase by adiabatic compression is canceled by 
efficient radiative cooling. Primordial molecules such as 
HD are particularly efficient for cooling at temperatures 
below 100 K. 

The prediction of the temperature runs in gravitational 
collapsing structures is highly non-trivial since the nega¬ 
tive specific heat of self-gravitating systems may lead to 
opposite behaviours than those in familiar thermal sys¬ 
tems. Observations of the galactic ISM show that unsta¬ 
ble gas tends to adopt a fractal structure over 4 — 5 orders 
of magnitude, and the smallest gas clouds may remain 
very c old and dense. In this context, Pfenniger & Combes 
( 1994 ) argue that cold H 2 gas may be a candidate for bary- 
onic dark matter provided that the smallest elements of 
the gas fractal structure, the clumpuscules, are sufficiently 
cold and compact. In their postulated conditions set es¬ 
sentially by the virial theorem, at low redshifts (T ss 3K, 
n « 10 ®cm“^) they noted that the conditions are met to 
form solid H 2 due to a much higher gas pressure than the 
one contemplated above in the unperturbed Hubble flow. 
However if the cold collapsed structures reach the CBR 
temperature from above, clearly condensed H 2 (in solid or 
liquid forms) cannot form below the redshift at which the 
radiation is at the H 2 critical point temperature, at 33.4 K. 
In other words, the largest redshift for the occurrence of 
condensed H 2 in equilibrium with the CBR in dense but 
cold clumps is at z = 33.4/2.726 — 1 = 11.25. 


5. Conclusions 

In the current understanding of the formation of the first 
bound structures during the dark age, the possibility that 
solid hydrogen flakes exist and modify the subsequent evo¬ 
lution must be considered. For regions in which the growth 
of inhomogeneities can be neglected, in widely different 
Universe models our chemical calculations show that be¬ 
tween redshifts 6 and 12 H 2 flakes may start to form. 
The subsequent chemistry and matter-radiation coupling 
should be significantly altered. Since several factors that 
can either damp or increase the formation of H 2 flakes, 
clearly one must take into account that the late dark age 
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thermochemistry is likely to be much more complex that 
earlier admitted. 

However at these redshifts most of the first struc¬ 
ture formation scenarios predict an already well advanced 
and widespread stage of reionization, which casts un¬ 
certainties about which processes win, or which frac¬ 
tion of the Universe between voids and collapsed regions 
are concerned by H 2 flake formation. The appearance of 
snowflakes at the end of the dark age would constitute the 
first “cosmic winter”. 
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Appendix 


We use the reactions rat es and Pij and photo-processes 
given by Galli & Palla (1998, 2002| ), except the hydrogen 
and deuterium recombination: 


H+ -b e" ^ H -b 7 

D+-be-^D-b7, (.20) 


with a reaction rate Orec (in cm s ^), which i s calc ulated 
from the approximation given by Abel et al. (1997): 


Q^rec 


exp 


-28.6130338-0.72411256 ln{T^) 


-2.02604473 X lO'^ 
-2.38086188 x 10“^ ln^(T„i) 
-3.21260521 X 10““^ ln'^(T„i) 
-1.42150291 X 10“^ ln'^(T„i) 
-b4.98910892 x 10“® ln®(Tm) 
-b5.75561414 x 10"^ In^(Tm) 
-1.85676704 x 10“® \n{Tm) 
-3.07113524 X 10“® ln®(T„i)j, 


(. 21 ) 


where is in eV. We consider also the reverse reaction 
(i.e., photoionization of neutral hydrogen and neutral deu¬ 
terium) : 

H-b 7 ^ H+ -be“ 

D-b 7 D+ -be- (.22) 


with the radiative rate coefficient (Jph (in s ^) through the 
CBR, defined by: 



diy 

exTp{hi'/ kTn) — 1 ’ 


(.23) 


where vth is the threshold energy for which photoioniza¬ 
tion is possible (here huth = 13.6 eV). We take the fre¬ 
quency dependent cross section tjph (in cm^) given by Abel 


et al. (1997): 


CTth 


6.30 X 10"^® 


(] 

4 exp 

4 — 4arctan(e)/e 



— exp(—27r/e) 


where e = 



(.24) 
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